Detection of XY behaviour in weakly anisotropic quantum antiferromagnets on the 

square lattice 
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We consider the Heisenberg antiferromagnet on the square lattice with 5 = 1/2 and very weak 
easy-plane exchange anisotropy; by means of the quantum Monte Carlo method, based on the 
continuous-time loop algorithm, we find that the thermodynamics of the model is highly sensitive 
to the presence of tiny anisotropics and is characterized by a crossover between isotropic and planar 
behaviour. We discuss the mechanism underlying the crossover phenomenon and show that it occurs 
at a temperature which is characteristic of the model. The expected Berezinskii-Kosterlitz-Thouless 
transition is observed below the crossover: a finite range of temperatures consequently opens for 
experimental detection of non-critical 2D XY behaviour. Direct comparison is made with uniform 
susceptibility data relative to the 5* = 1/2 layered antiferromagnet Sr2Cu02Cl2. 

PACS numbers: 75.10.Jm, 05.30.-d, 75.40.-s, 75.40.Cx 



Predictions of the Berezinskii-Kosterlitz-Thouless 
(BKT) theory [H for topological ordering with zero order 
parameter have been verified in many real systems, such 
as superfluid or superconducting films and Joseph- 
son junction arrays 0. However, despite the BKT the- 
ory being originally formulated as referred to 2D planar 
magnets, evidences of XY behaviour in real magnets are 
weak and limited to very peculiar cases . On the other 
hand, for S — 1/2 there exist several cuprous oxides that, 
besides being excellent realizations of 2D antiferromag- 
nets, are known to display a weak exchange easy-plane 
(EP) anisotropy: unambiguous observations of XY crit- 
ical behaviour in such systems are not easy to achieve, 
due to both the weakness of the anisotropy and the ex- 
istence of finite inter-layer coupling. As for the former 
point, the EP anisotropy observed in real magnets is 
usually 10""* ~ times the isotropic exchange cou- 

pling, and the signatures of BKT critical behaviour are 
often either too weak to be extracted from the isotropic 
thermodynamics or too close to the critical temperature 
to be experimentally accessible. The residual inter-layer 
coupling, even if orders of magnitudes smaller than the 
intra-layer one, drives the system towards a 3D transi- 
tion, which is actually triggered by the divergence of 2D 
intra-layer spin correlations. Therefore purely 2D critical 
behaviour of diverging quantities is most often masked by 
the onset of 3D long-range order. 

In this work we show that several non-diverging quan- 
tities are sensitive to the presence of EP anisotropy 
and display an evident and detectable crossover between 
isotropic and XY behaviour above the expected BKT 
transition. Such crossover occurs at a temperature which 
is characteristic of the model and is marked by peculiar 
features in the temperature dependence of non-critical 
observables. In particular, we present quantum Monte 



Carlo (QMC) data for the uniform susceptibility, the 
finite-size staggered out-of-plane magnetization, the spe- 
cific heat and the density of in-plane vortices. 

We consider the S' = 1/2 easy-plane antiferromagnet 
on the square lattice described by the Hamiltonian 
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where i — (ii, 22) runs over the sites of a square lattice, d 
connects each site to its four nearest neighbours, J > is 
the antiferromagnetic exchange coupling and A S (0, 1] 
is the EP anisotropy parameter. We use the reduced 
temperature t — T / J . 

The above model is studied by means of QMC sim- 
ulations based on the continuous-time loop algorithm 
H 013 for A = 0.001, 0.02, and 1, and on lattice sizes 
from L = 64 to L = 200. Each MC run consists of 10'* 
thcrmalization steps, followed by 1 ^ 1.5 x 10^ MC steps 
for measurements. Whenever possible, we implement im- 
proved estimators |^ for the quantities of interest. 

The class of EP antiferromagnets described by Eq. 
is characterized by the possibility for the model to sup- 
port non-linear topological excitations in the form of in- 
planc vortices (V) and antivortices (AV). Such excita- 
tions play a fundamental role in determining the ther- 
modynamics and critical behaviour of the system, as 
described by the BKT theory relative to the classical 
planar rotator model. When both quantum and out-of- 
plane fiuctuations are present, the BKT behaviour per- 
sists even in strongly quantum [S — 1/2) nearly isotropic 
(A ~ 10-3) models 8]. 

Nevertheless, there exist significant differences between 
the standard BKT phenomenology and the XY behaviour 
actually observed in quantum nearly isotropic antiferro- 
magnets, where the out-of-plane spin component plays 
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an essential role above the transition. 

Let us consider a strongly anisotropic EP model (A < 
1): At high temperatures the system, despite being dis- 
ordered, already supports topological excitations in the 
form of well separated V and AV in the easy plane. As 
t decreases V and AV attract each other, till V-AV pairs 
begin to form and a maximum in the specific heat is cor- 
respondingly observed. At the critical temperature i^j^^ 
V and AV are all paired and the topological transition 
occurs. 

In the nearly isotropic case (A <^ 1), the picture is 
modified as follows. At high temperature the system 
is isotropically disordered. When t decreases the EP 
anisotropy becomes effective enough to stabilize planar 
configurations and allow for V and AV to appear in 
the easy plane. This phenomenon, occurring at a tem- 
perature hereafter indicated as ico, may be thought of 
as a crossover between the isotropic and a more gen- 
uine XY behaviour, and suggests the temperature range 
'•BKT < ^ < ^co as the most appropriate for experimental 
observations. 

We consider the dimcnsionlcss uniform susceptibility 
X"", defined as 
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Fig. n shows its temperature dependence, which we find 
strongly characterized by the appearence of a minimum 
in the zz component. A similar minimum has been in- 
deed experimentally observed in Sr2Cu02Cl2 and sug- 
gested to be related to the onset of 2D XY behaviour [ll| . 

A sound argument in favour of this interpretation is the 
following. An infinitesimal uniform magnetic field causes 
a finite response in an antiferromagnet via the canting 
mechanism: adjacent spins antialigned in the plane per- 
pendicular to the field cant out of such plane and give rise 
to a net magnetization, while a negligible response comes 
from spins anti-aligned along the field direction. On this 
basis, we interpret Fig. At high temperature, Xu^ and 
behave in the same way and decrease upon decreas- 
ing t as antiferromagnetic coupling gets more effective 
and canting consequently harder. As t is further lowered, 
the anisotropy starts to play a role, and most of the spins 
antialign in the easy plane: therefore, compared to the 
isotropic case, the number of spins responding to a field 
applied in the plane (along z) decreases (increases). As a 
consequence, xS"^ decreases faster, while Xu^ slows down 
its decrease; the further reduction of out-of-plane fluctu- 
ations is eventually responsible for the low-temperature 
increase of Xu^- ^ clear minimum consequently appears 
and marks the crossover to XY behaviour at the tem- 
perature tco, where the out-of-plane component of the 
antiferromagnetic coupling becomes irrelevant. No par- 
ticular feature is instead seen at the transition. 

In order to check the direct relation between the 
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FIG. 1: Uniform susceptibility of the EP model for A = 0.02 
(triangles) and A = 0.001 (diamonds). Open (full) sym- 
bols represent the out-of-plane (in-plane) susceptibility. The 
stars represent the uniform susceptibility in the isotropic 
case 0,1^ Inset: finite-size magnetization of the EP 

model for A = 0.02 and L = 64 (diamonds), 96 (left trian- 
gles), and 128 (squares). Error bars are smaller than symbol 
sizes; arrows are the critical temperatures as estimated via 
finite-size scaling 0. 



crossover phenomenon and the out-of-plane fluctuations, 
we consider our data relative to finite-size out-of-plane 
staggered magnetization, l ; such quantity, which ob- 
viously vanishes at all temperatures in the thermody- 
namic limit, is quite useful when local spin configurations 
are under analysis, as in our case. In the inset of Fig. 
we show data for different lattice sizes and A = 0.02. By 
lowering t, increases (as expected in the isotropic 

behaviour) above tco, and decreases (as expected in XY 
behaviour) below tco- a stable maximum is seen at 
tco- From the above evidences we obtain the estimates 
tco(A = 0.02) = 0.30(1) and tco(A = 0.001) = 0.225(10). 
Signatures of the crossover are also present in the stag- 
gered out-of-plane susceptibility and out-of-plane corre- 
lation length (shown in Ref. 8), displaying a maximum 
at a temperature quite close to tco- 

The above results suggest that the crossover phe- 
nomenon is peculiar to the system and is due to the sup- 
pression of out-of-plane fluctuations. In particular, tco 
is the temperature below which the weakly anisotropic 
system behaves like a planar rotator model with a spin 
length effectively reduced by out-of-plane fluctuations. 
Other authors have referred to such crossover as due to a 
"spin-dimensionality reduction", meaning the loss of one 
spin component 
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By the semiclassical reasoning in Appendix B of Ref|3 
we can predict tco to depend on the anisotropy parameter 
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FIG. 2: Phase diagram of the weakly EP model on the square 
lattice. Down triangles: tco', up triangles: tsKT- Solid and 
dashed lines are logarithmic fits. Also reported are the tran- 
sition and crossover temperatures of Sr2Cu02Cl2 (see text). 
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where is the renormaUzed spin stiffness of the quantum 
isotropic model, and C is a constant. A logarithmic fit 
to our data, shown in Fig. |21 gives C = 160 and Pg = 
0.214 J, which compares well with the known value 
= 0.180 J. 

The onset of XY behaviour below ico is further sup- 
ported by the temperature dependence of the specific 
heat c. In the planar rotator model such quantity dis- 
plays a maximum well above the BKT transition. Fig. |2| 
shows our data for the weakly anisotropic model with 
A = 0.02: an embryonic peak, out of the isotropic com- 
ponent, is seen at t = 0.29(1), slightly below tco- To 
clarify the actual meaning of such peak, we have also 
considered the density of in-plane vortices with unitary 
vorticity, defined as 
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= -(i~8 s^s^ + 16 s^sp^sy^ 
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where m) defines a plaquette of the square lattice 

(indices are ordered counterclockwise). Choosing z as 
quantization axis, the above quantity is off-diagonal. The 
estimator of the bilinear term in Eq. can be found in 
Ref. In the same spirit we introduce the estimator for 
the quartic term in the context of the loop algorithm . 

By considering QMC data for A = 1, shown in the in- 
set of Fig. 121 we observe that the temperature derivative 
of p^ displays a maximum where the specific heat ex- 
hibits its peak. Both features mark therefore the onset of 
the formation of V-AV pairs. In Fig. Owe show that this 
picture is fully reproduced by our data for A = 0.02, thus 
confirming the XY character of our weakly anisotropic 
system. 

From the experimental point of view, the above find- 
ings may be firstly used to easily characterize the 



FIG. 3: Specific heat, vortex density and its temperature 
derivative for the EP model with A = 0.02. Dashed line: nu- 
merical derivative of splines for the vortex density data; solid 
line: specific heat of the isotropic model (from a numerical 
interpolation to the data of Ref. 0). Inset: same picture for 
A = 1 (vortex density data from Ref. IT3) . When not visible, 
the error bars are smaller than the symbol size; arrows are 
the BKT critical temperatures. 



anisotropy of the layered antiferromagnets: if a minimum 
in the out-of-plane component of the uniform susceptibil- 
ity is observed above the transition, this is a signature of 
EP anisotropy (while a minimum at the transition sug- 
gests an easy- axis anisotropy, as shown in Ref. 0). Once 
ico has been experimentally determined, Eq. ^ allows 
one to get an independent estimate of the bare anisotropy 
parameter A of Eq. (Q. As shown below, this is quite an 
essential point for interpreting the experimental data. 

We now consider the layered cuprate Sr2Cu02Cl2, 
whose intra-layer spin-spin coupling of the Cu^^ ions 
has been proposed to be governed by the Hamil- 
tonian with J = 1450 K. As for the anisotropy pa- 
rameter A, the experimental analysis hands us the 
renormalized value A'^'^p = 0.00014, extracted from low- 
temperature measurements of the spin gap in the out- 
of-plane branch of the 2D spin waves propagating in 
the Cu-layers, by means of the approximated expression 
G = 4JSZcV2A^, where Zc is the zero-temperature 
value of the spin-wave velocity renormalization coeffi- 
cient of the isotropic Heisenberg system Ig] . The above 
^oxp j^Q^ hence the bare value appearing in Eq. Q), as 
it already contains quantum renormalizations. In fact, 
a more refined self-consistent spin-wave-theory calcula- 
tion 01 leads to the expression G — AJSZc\/2Za^ with 
Za = 0.099, and hence A = A<=''p/Za =i 0.0014. 

Experimental data relative to the uniform susceptibil- 
ity of SraCuOzClz (from Refs. O and IT^ are shown 
in Fig. ^together with our results for A = 0.001; a 
constant offset for the experimental data has been intro- 
duced in order to take into account spurious temperature- 
independent contributions to the measured values of 
Xu"- The agreement is excellent for both the in-plane 
and the out-of-plane component; considering that no ad- 
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FIG. 4: Uniform susceptibility of Sr2Cu02Cl2 compared to 
the theoretical predictions for the EP model with A = 0.001. 
x's: Xu"^ of Sr2Cu02Cl2 (from Ref. [Tfl) (with gyromagnetic 
facto r =2); thin and thick +'s: xS^ of Sr2Cu02Cl2 (from 
Ref. llw with g—2 and g=2A6, respectively. Other symbols 
and error bars as in Fig. 

justment has been introduced for the temperature axis, 
the position of the minimum of is nicely reproduced, 
and a crossover temperature ico — 0.227 for the real 
compound is determined. When the critical region is ap- 
proached, experimental data deviate from the theoretical 
predictions due to the fact that a purely 2D model is no 
longer sufficient to capture the thermodynamic behaviour 
of the real magnet, which is in fact characterized by an 
experimentally observed Neel transition at = 0.176. 
Such transition temperature compares remarkably with 
the BKT transition temperature of the purely 2D model, 
'•BKT ~ 0.175(10): this confirms that the 3D ordering is 
induced by the incipient intra-layer BKT transition. 

From our analysis we conclude that the experimen- 
tal observation of a minimum in the transverse uniform 
susceptibility at a temperature which is well above the 
critical region, is a signature of EP character neatly 
detectable in real compounds even with very weak 
anisotropy. To this respect, we notice that the bare 
anisotropy value, needed for a quantitative compari- 
son between experimental and theoretical results, signif- 
icantly differs from the (measured) renormalized value, 
due to the relevance of quantum effects. Finally, two dis- 
tinct regimes, separated by the crossover, can be iden- 
tified in the disordered phase: an isotropic regime for 
t > tea and a disordered regime for <t <tco- It is 
relevant that the crossover temperature can be obtained 
by measuring non-diverging quantities above the critical 
region, which is an experimentally feasible task. 

In summary, by means of quantum Monte Carlo sim- 
ulations we have shown that two-dimensional S — 1/2 
square-lattice antiferromagnets with a weak easy-plane 
anisotropy display a crossover from a high-temperature 
isotropic behaviour to a genuinely 2D XY behaviour at 
a temperature tco which stays around 30% above the 



critical temperature tg^T ■ have moreover shown 

that Sr2Cu02Cl2 realizes these predictions, and there- 
fore stands as a very clean example of weakly planar 
quasi-2D antiferromagnet. This represent a clear evi- 
dence of non-critical XY behaviour in a real layered an- 
tiferromagnet. Other evidences for an isotropic-to-XY 
crossover in Sr2Cu02Cl2 come from neutron scattering 
[lll | and NMR experiments, and enforce our con- 

clusions. Further work on such crossover effect in the 
above compound, as well as in other layered antiferro- 
magnets with easy-plane anisotropy (as, e.g., Pr2Cu04), 
opens the perspective for a systematic investigation of 
the vortex phase in purely magnetic systems. 
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